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Qs ' The lateral and longitudinal profiles of hadronic showers detected by a prototype 

of the ATLAS Iron-Scintillator Tile Hadron Calorimeter have been investigated. 
This calorimeter uses a unique longitudinal configuration of scintillator tiles. Using 
a fine-grained pion beam scan at 100 GeV, a detailed picture of transverse shower 

O' behavior is obtained. The underlying radial energy densities for four depth seg- 

' ments and for the entire calorimeter have been reconstructed. A three-dimensional 

hadronic shower parametrization has been developed. The results presented here 
are useful for understanding the performance of iron-scintillator calorimeters, for 
developing fast simulations of hadronic showers, for many calorimetry problems 
requiring the integration of a shower energy deposition in a volume and for future 
calorimeter design. 



1 The Calorimeter 



0^ ' We report on an experimental study of hadronic shower profiles detected by 

^^ , the m-ototypc of the ATLAS Barrel Tile Hadron Calorimeter (Tile calorime- 

^ ' ter) El. The innovative design of this calorimeter, using longitudinal segmen- 

O , tation of active and passive layers provides an interesting system for the 

Oj' measurement of hadronic shower profiles. Specifically, we have studied the 

D , transverse development of hadronic showers using 100 GeV pion beams and 

i-Q ■ longitudinal development of hadronic showers using 20 - 300 GeV pion beams. 

The prototype Tile Calorimeter used for this study is composed of five mod- 
ules stacked in the Y direction, as shown in Fig. nl Each module spans 27r/64 
in the azimuthal angle, 100 cm in the Z direction, 180 cm in the X direc- 
j^ ■ tion (about 9 interaction lengths, A/, or about 80 effective radiation lengths, 

Xq), and has a front face of 100 x 20 cm^ □. The iron to scintillator ratio is 
4.67 : 1 by volume. The modules are divided into five segments along Z and 
they are also longitudinally segmented (along X) into four depth segments. 
The readout cells have a lateral dimensions of 200 mm along Z, and longi- 
tudinal dimensions of 300, 400, 500, 600 mm for depth segments 1-4. The 
calorimeter was placed on a scanning table that allowed movement in any di- 
rection. Upstream of the calorimeter, a trigger counter telescope (SI, S2, S3) 
was installed, defining a beam spot approximately 20 mm in diameter. Two 
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Figure 1. Schematic layout of the experimental setup (side view). 



delay-line wire chambers (BCl and BC2), each with {Z,Y) readout, allowed 
the impact point of beam particles on the calorimeter face to be reconstructed 
to better than ±1 mm u. "Muon walls" were placed behind (800x800 mm^), 
shown in Fig. |l| as "Back Muon Wall", and on the positive Z side (400x1150 
mm^), not seen in Fig. H^, of the calorimeter modules to measure longitudinal 
and lateral hadronic shower leakage. The data used for the study of lateral 
profiles were collected in 1995 during a special Z-scan run at the CERN SPS 
test beam. The calorimeter was exposed to 100 GeV negative pions at a 10° 
angle with varying impact points in the Z-range from —360 to 200 mm. A 
total of > 300,000 events have been analysed. The uniformity of the calorime- 
ter's response for this Z-scan is estimated to be 1% a. The data used for the 
study of longitudinal profiles were obtained using 20 - 300 GeV negative pions 
at a 20° angle and were also taken in 1995 during the same test beam run. 

2 Extracting the Underlying Radial Energy Density 



There are several methods for extracting the radial density $(r) from the 
measured distributions of energy depositions. One method was iised in the 
analysis of the data from the lead-scintillating fiber calorimeter U. Another 
method for extracting the radial density is to use the marginal density function 
f{z) which is related to the radial density <i>(r) u, 



*W = -Z3;j 



1 d 



f{z) dz^ 



^h^ 



We used the sum of three exponential functions to parameterize j(z) as 






(1) 



(2) 
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where z is the transverse coordinate, -Eq, a^, \ are free parameters, B = 
Si=i OiAi, X]i=i o-j = 1 and J_ f{z)dz = Eq. The radial density function, 
obtained by integration and differentiation of equation (|l|), is 

i— 1 

where Xq is the modified Bessel function. We define a column of five cells 
in a depth segment as a tower. Using the parametrization shown in equation 
(0) , we can show that the energy deposition in a tower can be written as 

E{z) ^Eo-^ Ell ««A, cosh(g) e-^, for \z\ < ^, (4) 

E{z) = f Eti a.A, sinh(2^) e"^, for \z\ > ^, (5) 

where h is the size of the front face of the tower along the z axis. 

3 Transverse Behaviour of Hadronic Showers 

Figure || shows the energy depositions in towers for depth segments 1-4 
as a function of the z coordinate of the center of the tower and for entire 
calorimeter. Here the coordinate system is linked to the incident particle 
direction where z = is the coordinate of the beam impact points at the 
calorimeter front face. With fine parallel displacements of the beam between 
—360 mm and 200 mm we expand the tower coordinate range from —760 mm 
to 600 mm. To avoid edge effects, we present tower energy depositions in the 
range from —650 mm to 500 mm. The tower energy depositions shown in 
Figures span a range of about three orders of magnitude. The plateau for 
\z\ < 100 mm {h/2) and the fall-off at large \z\ are apparent. We used the 
distributions in Figs. to extract the underlying marginal densities function 
for four depth segments of the calorimeter and for the entire calorimeter. The 
solid curves in these figures are the results of the fit with equations (4) and 
(™. The fits typically differ from the experimental distribution by less than 
5%. 

The parameters at and Xt, obtained by fitting, are listed in Table |lj. 
The tti and A, parameters demonstrate linear behaviour as a function of x: 
ai{x) = a.i + PiX and \i{x) = ^i -\- SiX. The values of the parameters a^. Pi, 
7i and Si are presented in Table ^ 
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Figure 2. Energy depositions of 100 GeV pions in towers of depth segments 1 — 4 as a 
function of the z coordinate: top left is for depth segment 1, top right is for depth segment 
2, middle left is for depth segment 3, middle right is for depth segment 4, bottom left is for 
over all calorimeter. Only statistical errors are shown. 



4 Radial Hadronic Shower Energy Density 

Using formula (0) and the values of the parameters a^, Ai, given in Table 0, 
we have determined the underlying radial hadronic shower energy density 
functions, <i>(r). The results are shown in Figured for depth segments 1-4 
and for the entire calorimeter. The contributions of the three terms of $(?•) 
are also shown. 

The function (f>(r) for the entire calorimeter has been compared with the 
one for the lead-scintillating fiber calorimeter of Ref. □, that has about the 
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Table 1. The parameters a^ and A^ obtained by fitting the transverse shower profiles for 
four depth segments and the entire calorimeter at 100 GeV. 



x{xin 


ai 


Ai(nini) 


02 


A2(nini) 


as 


A3 (mm) 


0.6 


0.88 ±0.07 


17 ±2 


0.12 ±0.07 


48 ±14 


0.004 ±0.002 


430 ± 240 


2.0 


0.79 ±0.06 


25 ±2 


0.20 ±0.06 


52 ±6 


0.014 ±0.006 


220 ± 40 


3.8 


0.69 ±0.03 


32 ±8 


0.28 ±0.03 


71 ±13 


0.029 ±0.005 


280 ± 30 


6.0 


0.41 ±0.05 


51 ±10 


0.52 ±0.06 


73 ±18 


0.07 ±0.03 


380 ± 140 


all 


0.78 ±0.08 


23 ±1 


0.20 ±0.08 


58 ±4 


0.015 ±0.004 


290 ± 40 



Table 2. The values of the parameters ai, Pi, 7^ and Si 





Oli 


P^ [l/K) 




7i (mm) 


5i (mm/A^) 


Oi 


0.99 ±0.06 


-0.088 ±0.015 


Ai 


13 ±2 


6±1 


02 


0.04 ±0.06 


0.071 ±0.015 


A2 


42 ±10 


6±4 


03 


-0.001 ±0.002 


0.008 ±0.002 


A3 


170 ± 80 


29 ±23 



same effective nuclear interaction length for pions u. The two radial den- 
sity functions are rather similar as seen in Fig. H (bottom right). The lead- 
scintillating fiber calorimeter density function <^{r) was obtained from a 80 
GeV 7r~ grid scan at an angle of 2° with respect to the fiber direction. For 
the sake of comparina_the radial density functions of the two calorimeters, 
the distribution from l3 was normalised to the <i>(r) of the Tile calorimeter. 
Precise agreement between these functions should not be expected because of 
the effect of the different absorber materials used in the two detectors, the 
values of e/h arc different, as is hadronic activity of showers because fewer 
neutrons are produced in iron than in lead. 

5 Radial Containment 



The parametrization of the radial density function, $(?"), was integrated 
to yield the shower containment as a function of the radius I{r) = Eq — 
-^ X]i=i ^i Ki{r/Xi), where Ki is the modified Bcssel function. The ap- 
proximation of the data for the radii of cylinders for given shower contain- 
ment (90%, 95% and 99%) as a function of depth with linear fits are r(90%) = 
(85±6)+(37±3)a;, r(95%) = (134±9) + (45±3)x, r(99%) = (349±7) + (22±2)a; 
(mm). The centers of depth segments, x, are given in units of A^*^. 
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Figure 3. Radial energy density, "J'(r), as a function of r for Tile calorimeter for depth 
segments 1 — 4: top left is for depth segment 1, top right is for depth segment 2, middle left is 
for depth segment 3, middle right is for depth segment 4, bottom left is for entire calorimeter. 
The solid lines are the energy densities <l?(r), the dashed lines are the contribution from 
the first term from, the dash-dotted lines are the contribution from the second term, the 
dotted lines are the contribution from the third term. Bottom right: Comparison of the 
radial energy densities as a function of r (in units of A^-*-* ) for Tile calorimeter (the solid 
line) and lead-scintillating fiber calorimeter (the dash-dotted line). 



6 Longitudinal Profile 

Our values of AE/Ax together with the data of 13 and Monte Carlo predictions 
(GEANT-FLUKA + MICAP) I are shown in Fig. | (top). The longitudinal 
energy deposition for our calorimeter is in good agreement with that of a 
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Figure 4. Top: The longitudinal profile (circles) of the hadronic shower at 100 GeV as a 
function of the longitudinal coordinate x in units of A^*^. Open triangles are data from the 
calorimeter of Ref.^ diamonds are the Monte Carlo (GEANT-FLUKA) predictions. The 
dash-dotted line is the fit by function (6), the solid line is calculated with function (7) with 
parameters from Ref.^ the dashed line is calculated with function (7) with parameters from 
Ref. The electromagnetic and hadronic components of the shower (crosses and squares), 
together with their fits using (6), are discussed in Section 8. Bottom: Longitudinal profiles 
of the hadronic showers from 20 GeV (open stars), 50 GeV (open squares) and 100 GeV 
(open circles) pions as a function of the longitudinal coordinate x in units of A/ for a 
conventional iron-scintillator calorimeter ''' and of 20 GeV (black stars), 50 GeV (black 
squares), 100 GeV (black circles), 150 GeV (black up triangles), 200 GeV (asterisks), 300 
GeV (black down triangles) for pions at 20° and of 100 GeV (black circles) for pions at 10° 
for Tile iron-scintillator calorimeter. The solid lines are calculated with function (7) with 
parameters from Ref.^ and the dashed lines are with parameters from Ref.* 



conventional iron-scintillator calorimeter. The longitudinal profile may be 
approximated using two parametrizations. The first form is 



dE{x) _ Ef (3- 



«+i 



-l3x 



dx T{a + 1) 



(6) 



where Ef = Eb^amr. and a and /3 are free parameters. Our data at 100 GeV 
and those of Ref. Q at 100 GeV were jointly fit to this expression; the fit is 
shown in Fig. |j (top). The second form is the analytical representation of the 



Y'KULCHITSKY: submitted to World Scientific on February 7, 2008 7 



longitudinal shower profile from the front of the calorimeter 



dx \ a \ Xq J \ \ Xj J Xq 

^ ' ' ^ e"^' li^i 1, a+1, (d-1) — , (7) 



a \^i J v' -^ 



V, 



where a, b, d, w are parameters, iFi is the confluent hypergeometric func- 
tion. Here the depth variable, x, is the depth in equivalent Fe, Xq is the 
radiation length in Fe and in this case A/ is A^*^. The normalisation factor 
iVp^ Ebeam/Xi T{a) {w Xq &~" + {1 — w) A/ d~°'). This form was suggested 
in 113 and derived by integration over the shower vertex positions of the longi- 
tudinal shower development from the shower origin. For the parametrization 
of longitudinal shower development from jthc shower origin, the well known 
parametrization suggested by Bock et al. H has been used. We compare the 
form (|^) to the experiraental points at 100 GeV using the parameters cal- 
culated in Refs. U and Q. Note that now we are not performing a fit but 
checking how well the general form (M) together with two sets of parameters 
for iron-scintillator calorimeters describe our data. As shown in Fig. H (top), 
both sets of parameters work rather well in describing the 100 GeV data. 
Turning next to the longitudinal shower development at different energies, 
in Fig. I (bottom) our values of AE/Ax for 20 - 300 GeV together with 
the data from Q are shown. The solid and xlashed lines are calculations with 
function (M) using parameters from 13 and a, respectively. Again, we observe 
reasonable agreement between our data and the corresponding data for con- 
ventional iron-scintillator calorimeter on one hand, andJbetween data and the 
parametrizations described above. Note that the fit in Q has been performed 
in the energy range from 10 to 140 GeV; hence the curves for 200 and 300 
GeV should be considered as extrapolations. It is not too surprising that at 
these energies the agreement is^ignificantly worse, particularly at 300 GeV. 
In contrast, the parameters of □ were derived from data spanning the range 
15 - 400 GeV, and are in much closer agreement with our data. 

7 The parametrization of Hadronic Showers 

The three-dimensional parametrization for spatial hadronic shower develop- 
ment is 

dE{x) Z/i=i TJ^ ^oixj^) 



dx 27r ^t,a,(a;)A,(x) 



Y'KULCHITSKY: submitted to World Scientific on February 7, 2008 8 



where dE{x)/dx, defined by equation Im), is the longitudinal energy deposi- 
tion. 

8 Electromagnetic Fraction of Hadronic Showers 




Figure 5. The fraction /^o in hadronic showers versus the beam energy. The star is 
our data, the sohd curve is the Groom parametrization, the dashed curve is the Wigmans 
parametrization, squares are the GEANT-CALOR predictions, circles are the GEANT- 
GHEISHA predictions and crosses are the CALOR predictions. 



Following cl, we assume that the electromagnetic part of a hadronic shower 
is the prominent central core, which in our case is the first term in the ex- 
pression (H) for the radial energy density function, <i>(r). Integrating f^^o over 
r we get /^o = aiAi/ X]i=i '^Ai- For the entire Tile calorimeter this value is 
(53±3)% at 100 GeV. The observed tt" fraction, f^^a, is related to the intrinsic 

e E' e/h-f'(E) 



e E'+h E' 



(e/h-l)-f' (E) + l^ 



actual fraction, /^o, by the equation f^o{E) 

where E'^^-^ and i?^ are the intrinsic electromagnetic and hadronic parts of 
shower energy, e and h are the coefficients of conversion of intrinsic electromag- 
netic and hadronic energies into observable signals, /^o = E'^mli^'em + -^D- 
jXhere are two analytic forms for the intrinsic^*' fraction suggested by Groom 
B I'^o{E) = 1 - (E/E^)"'-^ and Wigmans eI f'^o{E) = k ■ ln{E/E'^), where 
Eq — 1 GeV, m = 0.85 and k — 0.11_.We calculated /^o using the value 
e/h = 1.34 ± 0.03 for our calorimeter 13 and obtained the curves shown in 
Fig. 0. Our result at 100 GeV is compared to the modified Groom and Wig- 

aans parametrizations and to results from the Monte Carlo codes CALOR 
, GEANT-GEISHA and GEANT-CALOR (the latter code is an implemen- 
tation of CALOR89 differing from GEANT-FLUKA only for hadronic inter- 
actions below 10 GeV). As can be seen from Fig. H, our calculated value of 
f^o is about one standard deviation lower than two of the Monte Carlo re- 
sults and the Groom and Wigmans parametrizations. The fractions f^o (r) 
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for the entire calorimeter and for depth segments 1-3 amount to about 90% 
as r ^ and decrease to about 1% as r ^- X^^^ . However for depth segment 
4 the value of f^a (r) amounts to only 50% as r ^ and decreases slowly to 
about 10% as r -^ ^%^^ ■ The values of fTro[x) as a function of x are fitted by 
f^o{x) = (75 ± 2) - (8.4 ± 0.4)a; (%). Using the values of f^o{x) and energy 
depositions for various depth segments, we obtained the contributions from 
the electromagnetic and hadronic parts of hadronic showers in Fig. ^ (top). 
The curves represent a fit to the electromagnetic and hadronic components 
of the shower using equation (||). Ef is set equal to /^oi^beam for the elec- 
tromagnetic fraction and (1 — /^o)-E&eam for the hadronic fraction. The elec- 
tromagnetic component of a hadronic shower rise and decrease more rapidly 
than the hadronic one (agm = 1-4 ± 0.1, ah = 1.1 ± 0.1, Pem = 1-12 ± 0.04, 
Ph = 0.65 ± 0.05). The shower maximum position {x^ax = {ct/ P) A^*^) oc- 
curs at a shorter distance from the calorimeter front face (a;^^; = 1-23 X^^, 
1.85 A^*^). At depth segments greater than 4 A^^, the hadronic frac- 



x: 



h 



max 



tion of the shower begins to dominate. This is natural since the energy of the 
secondary hadrons is too low to permit significant pion production. 



9 Summary and Conclusions 

We have investigated the lateral development of hadronic showers using 100 
GeV pion beam data at an incidence angle of = 10° for impact points z in 
the range from —360 to 200 mm and the longitudinal development of hadronic 
showers using 20 - 300 GeV pion beams at an incidence angle of 8 = 20° . We 
have obtained for four depth segments and for the entire calorimeter: energy 
depositions in towers; underlying radial energy densities; the radii of cylinders 
for a given shower containment fraction; the fractions of the electromagnetic 
and hadronic parts of a shower; differential longitudinal energy deposition. 
The three-dimensional parametrization of hadronic showers that we obtained 
allows direct use in any application that requires volume integration of shower 
energy depositions and position reconstruction. 
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